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Abstract 

We present a new analysis of two-jet event shape distributions in soft collinear effective theory. 
Extending previous results, we observe that a large class of such distributions can be expressed 
in terms of vacuum matrix elements of operators in the effective theory. We match these matrix 
elements to the full theory in the two-jet limit without assuming factorization of the complete set 
of hadronic final states into independent sums over partonic collinear and soft states. We also 
briefly discuss the relationship of this approach to diagrammatic factorization in the full theory. 



* Electronic address: cwbauer@lbl.gov 

^ Electronic address: fleming@physics.arizona.edu 

^ Electronic address: clee@berkeley.edu 

§ Electronic address: sterman@insti.physics.sunysb.edu 

1 



1. INTRODUCTION 



Hadronic jets provide a window into the fundamental workings of quantum chromody- 
namics, since they contain within themselves the signatures of QCD at both weak and strong 
coupling. That hadronic final states in high-energy collisions are jet-like at all reflects the 
underlying parton- level perturbative interactions at weak coupling <^ 1, while the evolu- 
tion of individual partons into the final-state hadrons we detect depends on the transition 
to nonperturbative dynamics at strong coupling. Jets therefore test our understanding of 
QCD over a wide range of scales. 

At the same time that jets provide insights into QCD, they are key elements in signatures 
for new physics beyond the Standard Model. Unfortunately, the challenge of calculating 
QCD background events in the hadron collider environment of the Tevatron or the LHC 
remains formidable. To develop incrementally our ability to describe the structure and 
evolution of jets in these environments, it is useful to continue to develop the analysis of jets 
produced in the relatively cleaner environment of e'^e~ collisions. 

On the one hand, one can study observables that depend on specifying the actual number 
of jets in the final state by defining a jet algorithm, which sets criteria for what constitutes 
a jet. On the other hand, one can extract much useful information about the structure of 
the final state from simpler observables called event shapes, which do not depend on a jet 
algorithm, but are simple functions of the momenta of all the particles in the final state. 
Two-jet event shapes e are those whose distribution near e = is dominated by events with 
two nearly back-to-back collimated jets of particles. This is the case for e+e~ annihilation at 
high energy, and there is no need to specify precisely what constitutes a jet. Some familiar 
examples are thrust r = 1 — T Hj, jet broadening B [3|, jet masses 0, 0, 0], and the 
C-parameter A generic class of event shapes for events e~^e~ ^ X at center-of-mass 
energy Q may be defined by 

where rji and are the rapidity and transverse momentum of final state particle i with 
respect to the thrust axis of the event. The function fe is sufficiently well behaved to 
guarantee infrared safety. Choices of fe that give familiar event shapes are 

/.(r/)=e-l''l, fB{v) = l, fc{v) = ^- (2) 

coshr/ 

A recently-introduced class of event shapes r^, the angularities [sl, 0], for which 

/a(r])=e-l''l(i-'^), (3) 

generalizes thrust and jet broadening. The thrust corresponds to a = 0, and the broadening 
to a = 1. Any a with — oo < a < 2 defines an infrared-safe observable. Studying how 
the behavior of the angularity distributions vary with a provides insight into the intrinsic 
structure of jets jl, 1^ . 



In the endpoint region of two-jet event shape distributions, fixed-order perturbation the- 
ory alone is insufficient to make accurate predictions. This region is dominated by pertur- 
bative logarithms of e and by nonperturbative power corrections. A tool to separate the 
perturbative and nonperturbative contributions and set up a resummation of the logarithms 
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in these distributions is factorization, which separates the effects of the hard scattering, jet 
evolution, and hadronization occurring at different length and energy scales. An observable 
factorizes if it can be calculated as the convolution of perturbatively-calculable functions 
and nonperturbative functions, which are typically universal in the sense that the same non- 
perturbative function contributes to multiple processes. A two-jet event shape distribution 
that factorizes takes the form 

1 da 

(To de 



H{fi) j deide2 de^ 5(e - ei - 62 - e^) Ji(ei; /i) J2(e2; yu)5'(es; /i). 



(4) 



where H{fi) is a perturbatively-calculable hard function, and Jj and S are jet and soft 
functions respectively. This form summarizes the leading-power behavior, 1/e x logarithms 
of e, of the cross-section. Corrections are suppressed by powers of e. 

It is important to emphasize that the precise definitions of the jet and soft functions 
are somewhat arbitrary, and differ from treatment to treatment. Analyzing both of these 
functions perturbatively, one finds logarithmic dependence on the event shape variables Cj 
and Cs- This dependence can be resummed, from which we can generate the resummed 
dependence of the cross-section on e. The original resummations of thrust and related event 
shape variables [ll[ used the feature that the entire next-to-leading logarithmic behavior is 
naturally absorbed into the jet functions. These discussions do not need to include a soft 
function. The general structure of Eq. (jlj), including both soft and jet functions, and its 
relationship to Sudakov resummation at arbitrary logarithmic level was explored in a general 
context in Ref. T2|- Perturbative applications to jet shapes were described in more detail, 
with explicit constructions for the jet and soft functions in Ref. j^. These treatments assume 
that perturbation theory provides a leading-power description of infrared-safe observables in 
the manner of an asymptotic series, with power corrections whose behavior can be inferred 
from the ^structure of the perturbative series at high orders, and/or the running of the 



coupling [13j. Both the jet and soft functions also generally receive power corrections in 
this manner, and nonperturbative corrections to event shapes based on_both fixed-order and 
resummed cross sections have been widely discussed, and reviewed in 
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Recently the factorization of jet cross sections and event shape distributions near the two- 
jet kinematic eiidpoint has been revisited in the language of soft coUinear effective th eory 



(SCET) [15|,ll6|,ll7|,ll8 



3 Hill, [22 



This approach provides an elegant way, developed in 
to reproduce the results of the traditional QCD factorization at leading power, and provides 

' 3, I2I, 2^ and nonperturbative 



a framework to organize perturbative resummations [23 



power corrections [21 



case of massive jets [28 



27 
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These analyses of massless jets have also been extended to the 
. The effective field theory also provides a systematic framework 
to move beyond leading-power results. In SCET, the factorization in Eq. (jlj) follows from the 
usual separation of short- and long-distance physics in effective field theories, as well as from 
the decoupling of interactions between collinear and soft fields in the leading-order SCET 
Lagrangian. The hard function H is the square of a Wilson coefficient from the matching 
between QCD currents and SCET operators, the jet functions Jj are matrix elements of 
operators built of collinear fields, and the soft function S is the matrix element of an operator 
built from ultrasoft (usoft) fields. 

One goal of the SCET formulations of factorization in [l^, H, [2l|, H^, 28 1 has been to 
use effective theory methods to improve our insight into nonperturbative corrections by 
analyzing matrix elements directly in terms of hadronic states. As we shall see, the jet cross 
sections that we study below can be reformulated in terms of hadronic matrix elements, 
involving nonlocal products of currents and the energy- momentum tensor 30|, l3l|]- There is 
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an important subtlety, however, in the apphcation of effective field theory methods to semi- 
inclusive sums over asymptotic states. This is the assumption that hadronic final states X 
can be factorized consistently as 

|X) = |X,) iXsoft) (5) 

where Xc and Xsoft are collinear and soft final states, which carry non-singlet color in gen- 
eral. Indeed, it is inconsistent to impose a color-singlet condition on the collinear and soft 
states individually [13]. Furthermore, the relation Eq. (j5]) must be associated with power 
corrections that in principle depend on the observable to which we intend to apply it. Thus, 
corrections to Eq. are not well-defined without additional input. Similar limitations 
apply to an appeal to parton-hadron duality to justify Eq. (jS]), since this principle does not 
come with a systematic method of estimating corrections. Finally, if we were to take both 
sides of the relation literally, it would have to hold at infinite times, which would violate the 
spirit of matching the effective to full theory at a perturbative scale, and only then evolving 
to nonperturbative scales. These considerations motivate our analysis below. 

In this paper, we show that the assumption Eq. is not necessary to prove the factor- 
ization Eq. (jlj). We show that the phase space delta function 5(e — e(X)), which restricts the 
final states in da/de to the events with e(X) = e, can be expressed as an operator S{e — e) 
acting on the final state X. The operator e, which can be constructed in quantum field 
theory, depends on an energy flow operator as well as an operator t, which picks out the 
thrust axis of a final state X. Using these ingredients to express S{e — e(X)) as an operator, 
we remove its dependence on the final state X and so are able to sum over the complete set 
of states before it ever becomes necessary to assume the factorization of states, Eq. ([5]). 

In the next section, we demonstrate the factorization of the event shape distribution 
da/de in SCET, making use of the event shape operator e. In Section [3] we construct the 
energy flow operator that enters in the construction of e and discuss the determination 
within SCET of the thrust axis on which e depends. Section H] contains our conclusions. 



2. FACTORIZATION PROOF 

A. Event shape distributions in QCD 

We start with the distribution in the event shape e in e~^e~ — *• hadrons, which is given in 
full QCD by 

^ = ^ E ^ ^) r (2^)'^'(9 - PxMe - e{X)) , (6) 

where q is the total incoming momentum, with q"^ = Q"^. Here X labels the hadronic final 
state, and e(X) denotes the value of the event shape variable for a given final state X. 
To leading order in the electroweak couplings, the process e~^e~ ^ X is mediated by the 
partonic s-channel transition e~^e~ — > qq, with an intermediate photon or Z boson. The 
leptonic part of this partonic transition matrix element can be calculated, and one finds 

\M{e^e- ^X)\'=J2 (0| 1^) (^1 jT |0) , (7) 

i=V,A 
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where we have defined the vector and axial currents, 



Ji — Hf^ i If ■> 

with Ty = 7^ and = 7^7^. The leptonic tensor is given by 



-V 



3g2 



Q2 + 



/ 



(Q2 _ M|)2 



(8) 

(9a) 

(9b) 



where fermion / has electric charge in units of e, and vector and axial charges Vf,af 
given by 



1 



2 sin 9]y cos 9w 



{Tf - 2Qf sin^ 9 



w 



Tf. 

2 sin 9w cos 9w ■' 



(10) 



In Eq. ([H]) a sum over colors a and fiavors / is understood. 

Writing the four-momentum conserving delta function as the integral of an exponential, 
and using the dependence on px in the exponential to translate one of the two currents to 
the position x we can write the distribution as 



da 

de 



^,Y.I ^'^^^''^ E ^;.(0|jf (^) W(^|jr(0)|0)5(e-e(X)), (11) 
^ X i=V,A 



B. Eliminating the dependence on the final state 

The delta function 5{e — e{X)) restricts the sum over final states to those states giving 
the same value e of the observable event shape. This means that we cannot perform the sum 
over the complete set of final states. However, as we will now show, it is possible to write 
the event shape e(X) as the eigenvalue of an operator acting on the final state X. This 
can be achieved using the definition of the transverse energy fiow operator St{ti), which was 
introduced in |21|]. Its action on a hadronic state X is given by 

ST{v)\X) = J2\pJ\^iv-V^)\X), (12) 

where pf is the transverse momentum of the ith particle with respect to the thrust axis, and 
Tji is the rapidity of the ith particle. The thrust axis is defined to be the unit vector t which 
maximizes the sum ^j|pj-t|. In the event shapes of Eq. ([T]), rapidities and transverse 
momenta are measured with respect to this axis. Thus implicit in the action of Exiv) 
\X) is the determination of this thrust axis t(X). Using the energy flow operator we define 
an operator e, which returns the value of the event shape for a given state X, 

1 1"^ 

e\X)^e{X)\X) = - j d7^av)^T{v;t)\X) , (13) 

where t is an operator that returns the value of the thrust axis t{X) when acting on the final 
state X, and we have denoted exphcitly the dependence of Srirj) on this axis. Although the 
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present argument does not rely on explicitly constructing t, it is nevertheless possible to do 
so, as we show in Sec. 13 B 1[ In Sec. 13 B 21 we argue that in SCET we can choose the thrust 
axis to be in the jet direction n appearing in the two-jet current, so that no t operator need 
act on the final state at all. 

Using the thrust axis and event shape operators t and e, we can remove all dependence 
on the final state in the factor 6{e — e{X)) in Eq. ffTTl) and can therefore perform the sum 
over the complete set of final states. This gives 



da 1 



de 2Q2 



[ d'xe^^- Yl (Olif (^)^(e - ^)jr(O) |0) ■ (14) 

i=V,A 



The expression above involves a delta function of the operator e, which requires further 
comment. Heuristically this delta function is a way of treating the factorization of all 
moments of e at the same time. To see this we first define the delta function operator 
S{e — e) through a Taylor series expansion in e: 



S{e - e) = 6{e) + e6^^\e) + ■ ■ ■ + ^ 6^"\e) + ■ ■ ■ . (15) 



From this expression it is clear that the nth term in the series is the nth moment of the event 
shape distribution. Thus if we integrate Eq. f|T^ against e" the delta function operator on 
the right picks out the nth moment of the event shape distribution. 

In order to factorize this matrix element, we need to match the full theory currents onto 
operators in SCET, and to construct explicitly the operator e in SCET. The operator St{i]) 



is related to the energy-momentum tensor by [30, l31 



St{7]) = ^ / d0 lim / dtniToi{t,Rn) . (16) 

cosh 7] Jo Jo 

In Sec. 13 A\ we will prove Eq. (ITB]) using the energy- momentum tensor T^,^ written in terms 
of fields corresponding to the hadrons in the state X. In the proof of factorization below, we 
will instead use its presumably equivalent form in terms of quark and gluon fields in QCD 
and SCET. We are free to use either form as an operator is independent of its representation. 



C. Matching onto SCET 

We now match the currents and energy fiow operator onto SCET. The discussion in this 
section is purposely kept brief, and for more details of the techniques used and for definitions 



of our notation we refer the reader to Refs. [20|, l2l|, |28|. To reproduce the endpoint region of 
the two-jet event shape distribution, we match the QCD currents jf onto SCET operators 
containing fields in only two collinear directions: 

Jiix) = Y YCn,n2{PuP2;fJ')On,n2{x;PuP2) ■ (17) 
ni,n2 pi,P2 

The operator O can depend on the label directions rii and n2, as well as the label momenta 
Pi,P2- Recall that, in SCET, collinear momenta p'^ = p'^ + are divided into a large label 
piece, = (n ■p)n^/2 and a residual piece, k^, where n -p is 0{Q), p± is 0{QX), and 
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k is (9(QA^). The residual momenta are the same size as ultrasoft momenta, fcus, of 0{QX^). 
The small parameter A is of order \/Aqcd/Q- After the Bauer-Pirjol-Stewart (BPS) field 
redefinition with usoft Wilson lines the current in SCET is given by 

^711712 (-^i Pi) P2) 6 ''^ ^ Xni,pi{x)Yni{x)'^ iY n2{^)Xn2,P2i-^) 1 i^^) 

where Ty = 7^ and = 7^75. In Eq. (1171) we sum over directions rij of the light-cone 
vectors Ui = (l,nj). In the center-of-mass frame the two jet directions are constrained to 
be back-to-back, which eliminates the sum over one of the directions rij, leaving a sum over 
a single n. The field redefinitions replace all usoft-coUinear interactions with usoft Wilson 
lines, Yn^,Yni, in the directions and color representations of the corresponding collinear 
fields. 

Since collinear fields in different directions do not couple to one another in SCET at 
leading order, the forward matrix element in Eq. f|T^ vanishes, unless the directions n of 
the operator O and its complex conjugate agree with one another. By the same argument, 
the label momenta on these fields and their complex conjugates have to agree with one 
another. This gives 

% E C'nn(Pl,P2;/i)C„n(Pl,P2; 

n Pi,P2 

i=V,A 

where we have made all color indices explicit. Finally, we demand that the label momenta 
that appear in the exponentials of this relation equal the total momentum, thus requiring 
n ■ pi = —n ■ p2 = Q and Pi = P2 = 0. 



/i) j d^x eiQ^°ei(P2-Pi)-x 



{x)6{e - e) kip Y^^^KYnXi J (0) |0) 



(19) 



D. SCET scaling and the event shapes 

For our reasoning below, it will be important to estimate the contributions to event 
shapes from the label and residual momenta of collinear particles in the final state as well 
as from usoft particles. To be concrete, we consider the angularities e = r^, from Eq. ([1]), 
for which fa{v) = e"'^"-^"'^-*, and follow the logic of jsj. The contribution of an individual 
particle to Ta{X) is 

l-a 



mm 



For usoft particles we have px ~ QA^ and the ratio in parentheses i? ~ 1, while for collinear 
particles pt ~ QA and R ~ A^. The nominal contribution of a usoft particle to is thus A^, 
while for a collinear particle it is A^~". This is the case for any value of parameter A -C 1, 
For a < 0, the scaling behavior of the event shapes is A^, independent of a, and dominated 
by the momenta of usoft particles, independent of the collinear particles. For < a < 1, 
the event shapes scale as A^~°. We will establish our results below for all angularities with 
a < 1 and related event shapes, keeping all contributions to the event shape which are at 
least of order A^. 

For a > 1, the contribution of "ultra-coUinear" particles, for which pt ~ QX^ and i? ~ A^ 



become important 2l|. Their contribution to the event shape is of the order A^*^^ "\ which 
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is at least as large as for a > 1. Such particles are described by collinear modes in SCET 
with purely residual (zero label) transverse momenta. Our analysis below applies to jets 
whose typical constituents have nonzero transverse label momenta, and so is appropriate 
only for a < 1. 

The subdominance of ultra-coUinear contributions for a < 1 corresponds to the observa- 
tion that at a given px, wide-angle emission with energy comparable to pt contributes to 
the thrust and related event shapes more than does collinear emission with energies much 
larger than pt 32| at the level of both perturbative and nonperturbative corrections. 



E. Factorizing the vacuum matrix element 

At this point, the field content in the forward matrix element is completely factorized. 
To complete the factorization of the matrix element in Eq. flTI?]) . we must in addition show 
that the operator 6{e — e) also factorizes into collinear and soft parts. From Eqs. f[T^ and 
flTB]) we can see that the operator e can be written in terms of the energy-momentum tensor, 
as well as the operator t, which determines the thrust axis. As we show in Sec. 13 Bt the 
thrust axis is determined solely by the labels on the fields in the two-jet operator Onn, for 
which one simply finds that t = n. We will also use our observations in Sec. 12 Dl above, that 
the event shape values, and hence the action of the operator e on final state, are dominated 
by the contributions of usoft momenta, collinear label momenta fii-p and p±, and only one 
component Ui-k of collinear residual momenta, for particles in each collinear direction nj. 

The energy-momentum tensor, which is defined as 

where the sum is over fields 0j in the Lagrangian C, simplifies in SCET since after the BPS 
field redefinition the Lagrangian separates into 

As a result the energy-momentum tensor is a direct sum over contributions from fields in 
the different sectors. We must remember, of course, that Eq. fl2^ holds only at leading 
order in A in SCET, and that there are power-suppressed terms in the SCET Lagrangian 
in which interactions between collinear and usoft fields do not decouple following the BPS 
field redefintion 33, 34, 35, 36l . 



Using the definition of the event shape operator given in Eq. (fT3|) . we find the important 
result in SCET 

e = e^, -|- 6fi, -|- , (23) 

where 

1 /"^ 

^'^q/ ^^/e(^)^r(^), (24) 

and S^{ri) is defined using Eq. (fT6|) . but using the energy-momentum tensor derived from 
£j. This means that we can write the delta function constraining the value of the event 
shape to its observed value as 

5{e- e) = J (ie„5(e„ - e„) J dcnS^en - en) J deusSie^s - Cus) (5(e - e„ - - e^s) ■ (25) 
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Finally, we use the fact that the operators (e^) are constructed only from collinear fields 
in the n (n) direction and Cus only from usoft fields. Thus 

[en, Xn] = [Cfi, Xn] = [e„, Y] = [en, Y] = [Cus, Xn] = [ens, Xn]=0. (26) 

This enables us to rewrite Eq. (fT9|) in factorized form, separating the vacuum expectation 
values of mutually commuting fields, 

dcT 1 ft 

-j^ = ^Q^^\^nn{Q,~Q; J d'^x J dendendesS{e- en - en - Cs) (27) 

X ^ Tr (0| Xn,Q{x)i3S{en - e„)Xn,Q(0)^ |0) Tr (0| Xn-Q{x)J{en - es)Xn -q(0)5 |0) 
xTr(0|Fl(a;)y;t(a;)5(e,-e.)Fn(0)Fn(0)|0) L\TfUiT,,),s , 

i=V,A 

where L* = g^'^U^^, the traces are over colors, and we now make the spin indices explicit. 
The collinear matrix elements define jet functions J'n^n according to 

— Tr (0| Xn,Q{x),5{en - e„)Xn,Q(0), |0) ^ j —^^(^^^^"''^ Men, k^; ^) (^f| (28a) 

-l-Tr{0\xn,-Qix)Jien-en)xn,^QiO)5\0) = J ^^^^^e"^''^ ^^(en, T ; /i) , (28b) 
while the usoft matrix element defines a soft function 

^ Tr (0| Yn{x)Y^{x)S{es - e,)y„(0)F,(0) |0) ^ J ^e--^(e„ r; f,) . (29) 

Furthermore, we can use that J'n{en, fi) depends only on single light-cone component 
= n-k oi the residual momentum, as only in ■ d appears in the n-coUinear SCET La- 



grangian [l7j at leading-order in A. Similarly J7s(e.ra, /~; /i) depends only on l~ = n-l. To 
express the cross section, Eq. (1271) . in terms of the jet and soft functions directly, we follow 
a variation of the reasoning of Ref. j28| for massive quarks. 

Given the discussion of Sec. 12 Dl above, the residual transverse {k± and /_l) and parallel- 
moving {k~ and components of the jet functions enter the angularity event shapes only 
at nonleading power in A, as long as the parameter a < 1. This is the case for all values of 
A ^ 1, and therefore holds for perturbative, logarithmic corrections to the angularity cross 
sections, as well as in the nonperturbative region. For the following discussion, we therefore 
define the jet event shape operators e„, and en to measure only the residual components k~^ 
and /~ in the collinear sector. Corrections to this approximation from extending integrations 
over residual momenta to order Q are perturbative and nonlogarithmic, and hence can be 
absorbed into the matching coefficients Cnn- For simplicity, however, we do not change our 
notation for the matching coefficients to reflect this here. 

Because the jet functions are independent of the residual transverse momenta, it is natural 
to change variables to their sum and difference: K± = k±+l±, k± = {k±—l±)/2. The integral 
over the sum gives {27i)'^6'^{xj_). In this notation, the formula Eq. fl271) for the event shape 
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distribution can now be written as 

^^L^ J £kx_ \CnniQ,-Q;i^)\'^ J d^x j dendendes6{e- en - Cn- Cs 

dk^dl' 



de 



X 6 



X 



6 



X 



4(27r)4 J (277) 



d r 



-i{r++k+)x~ /2-i{r~+l-)x+/2-irj_-Xj 



where the factor L is defined 



L = L^Tt 



+ ^^Trf^7^5k^5 



(30) 



(31) 



In Eq. ( l30l) we have integrated over k',l^, and K± to generate deha functions setting all 
components of x to zero. This allows us to perform the integrals over x. 

In Eq. ( 1301) there remains a sum over label directions n and an integral over the residual 



momenta k±, which combined are simply an integral over total solid angle [28 



4 



(32) 



where the overall factor arises from the magnitude of the label three-momentum of the jet 
in direction n, which is |p| = Q/2. We carry out the integral over solid angle to obtain 
our final result. To simplify our expression we define the Wilson coefficient C2{Q'-, ^) = 
CnniQ, —Q', f^) since it is independent of n. We now define jet and soft functions integrated 
over all momenta: 



Jnifini A'') — 
Jnifini A'') = 

S{es;fi) = 



'^Jni^^ni k 1 /i) 



dl 1- \ 

Zn 

d^T 

S'(es,r;/i) . 



(277)' 



(33a) 
(33b) 
(33c) 



We note that the differential cross section can be expressed in terms of the total Born 
cross-section for e^e" qq 



^0 = 2^ 



3g2 



/ 



2\1 



{Q 



2 - M|)2 



(34) 



and we find as the final result for the differential event shape distribution: 
1 da 



(To de 



|C,(Q; ,)f j ie„ de, ie. S(e - e„ - e, - e.) J„(e„; ,)Me,^ ,)S(e.. ,) 



(35) 



This is the main result of this paper, in which the hard scattering function H{^) of Eq. (jl]) 
is identified with the absolute square of an SCET matching coefficient, and the jet and 
soft functions have been given effective theory definitions. Corrections suppressed at least 
by a power of A to this formula enter from effects at subleading order in the SCET power 
expansion. 
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F. Comment on final states in SCET 



In earlier studies of event shape distributions and two-jet cross sections in SCET 19|, l20|, 
21 . 25 . 28| . the hadronic final states X appearing in the cross-sections were not summed 
over before matching the full and effective theories. In these cases, the state X was itself 
factored into separate coUinear and soft parts 

\X) = \Xn) \X^s) . (36) 

In these studies the factorization theorem Eq. fl5Sl) is of the same form, but the jet function 
J'n is defined by 



(0|x„,Q(x)^|X„)(X„|x„,Q(0)^|0)5(e„-e(X„)) = / ^e'^'^-^ J-„(e„, /i) 



(37) 

and similarly for and S. Such definitions, of course, can be derived from the vacuum 
matrix elements in Eqs. ( 128|) and ( l29i) . if we assume that each set of XniXn^X^s form 
individually a complete set of states. We note, however, that the fields Xn,n overlap with 
states Xn,n only if these states are color triplet. Similarly, the states X^ in general carry 
non-singlet color. At this level, discussions that begin with the assumption of Eq. (136|) 
yield results that are equivalent to those of our reasoning. We have outlined the conceptual 
drawbacks of the use of this assumption in the Introduction. The advantages of avoiding 
the separation of states as in Eq. (!36|) . however, are not only conceptual. By defining the jet 
functions in terms of matrix elements as in Eqs. fl28l) and fl29l) . we can estimate corrections to 
the factorized cross section. These are due on the one hand to the leading-A approximation 
in the equation for the QCD Lagrange density, fl22]) . and on the other hand to our ability to 



calculate or otherwise estimate through power corrections j2l|] the specific vacuum matrix 
elements involving the operators Cj. 

We may also compare the role of partonic final states in SCET to their treatment in 
Ref. Wi, in which factorized iet and soft functions are each associated with individual sums 



over final states. In Ref. |8|, as in previous derivations of resummed event shapes [3 



381 ]. final states are always partonic. The factorization is carried out at the level of the 
diagrammatic cross section order-by-order in the coupling and region-by-region in phase 
space. Each such region is characterized by a definite set of jet and soft final-state partons, 
which are grouped into factorized jet and soft functions. In Ref. jsl, each of these factors is 
initially associated with a limited phase space for its final-state partons. These initial sums 
over partonic final states are then extended into unconstrained phase space sums for the jet 
and final states by redefining the functions systematically to avoid double counting in any 
region that is infrared sensitive. This procedure relies on the exponentiation properties of 
eikonal cross sections fs^, 0, 41. 4^ . 



In the diagrammatic treatment just described, specific classes of nonperturbative power 
corrections may be inferred from ambiguities in the resummed expressions fl3|]. A full 
treatment of nonperturbative corrections would require analyzing the observables in terms 



of matrix elements of the full theory, as advocated in [21|, |30|, |31|, a viewpoint that is 
consistent with the SCET analysis described above. 
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3. DEFINITION OF REQUIRED OPERATORS 



In this section we explicitly construct the operator e defined in Eq. f[T3|) . which we used to 
eliminate the dependence on the final hadronic state X of the weight 5{e — e(X)) appearing 
in the event shape distribution da/de. This operator e itself depends on the transverse 
momentum fiow operator ST{r)) and the thrust axis operator t. We first confirm Eq. flTB]) 
giving St{ti) in terms of the energy- momentum tensor T^,^, and then argue that t can be 
replaced with an axis that depends only on labels in the effective theory two-jet operator 
up to subleading power corrections (at least for a large set of event shapes which we identify). 

A. Energy Flow from Energy-Momentum Tensor 

The transverse momentum flow operator Sxirj, t) in Eq. f[T51) is related to the energy flow 



operator S{n) deflned in Ref. [43|, whose action on states |X) = \ki, . . . , kn) is 



S{h)\X)=J2^^S'in-h,)\X), (38) 

where uji is the energy of particle i, n is the unit vector pointing in the direction [6, 0), and 
6'^{h — hi) = 6{cos6 — cos6i)6{(j) — (pi)- To change variables from 6 to r], we use the relation 
cos 6 = tanh rj for massless particles, to obtain 

5(cos 9 — cos 9i) = 5(tanh r] — tanh r]i) = cosh^ r] 5{r] — rji) . (39) 

We deflne 9 and 7] to be measured with respect to the thrust axis. Also, the energy uji and 
the momentum |p^| transverse to the thrust axis are related by 

Ipf I = cui sin 9i = . (40) 

cosh rji 



Thus we can relate the transverse momentum flow operator in Eq. f[T2|) to the energy flow 
operator in Eq. (1551) by 



]^ r2TT 



^t{v) = —rr- / d<pS{h). (41) 
cosh f] Jo 

We will now verify that S{h) is related to the energy- momentum tensor by 

POO 

£{h) = lim / dt hiToi{t,Rh) , (42) 



the form of which was introduced in Ref. [3l|.^ We will think of the variable i? as a truly 
macroscopic distance, literally the distance from the scattering interaction to a detector. 
Thus, R is typically many orders of magnitude larger than the typical inverse minimum 
mass scale of the process. 



^ A similar formula integrating over distance in the direction fi with time taken to infinity is derived in 
Refs. ^\m\- 
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Consider the observable 

F{h,q) = J2 [ d'xe''^-^ {0\j\x)S{n) \X) {X\j{0) |0) , (43) 

where the sum is over a complete set of out states X, and the current j{x) couples to the 
field 0. Using the definition of the energy flow operators given in Eq. fl38|) . one can easily 
show that this observable gives the weighted cross section 

Fin,q)= (27r)^5^(px-g) 5^c.,52(n-n,)(0|jt(0)|X)(X|j(0)|0) , (44) 

states X particles i(^X 

where is a unit vector in the direction of the three- momentum of particle i. Our aim is to 
show that the representation of the energy flow operator in terms of the energy-momentum 
tensor, as given in Eq. P2|) . reproduces this result, with corrections entering with inverse 
powers of the distance R from the interaction region to the detector. 

To begin, we insert another sum over out states X' between the current j{x) and S in 
Eq. (113]). This gives 

states X,X' 

We then observe that since S{n) is at the macroscopic distance R from the scattering, all 
hadrons in the states X, X' will have stopped interacting by the time they reach the position 
of the operator S{n). The matrix elements of Sl-h) between states X, X' are thus those of a 
free hadronic theory. 

To prove Eq. fH2|) . we must show that S{n) acts on these hadronic states according to 
Eq. (!38|) . We do so by plugging the energy-momentum tensor T^^ appropriate for hadrons of 
a given type into Eq. fH2]) and testing its action on these states in the appropriate free fleld 
theory. This works for hadrons of any spin. Below we will demonstrate this explicitly for 
real scalars and Dirac fermions. Note that the scalar and Dirac flelds represent the hadronic 
flnal states, not the partonic states. Since the hadrons are non-interacting, we only need to 
consider free field theories. 



1 . Scalar fields 



We will first evaluate the matrix elements of S{n) in Eq. fH5l) for a free, neutral scalar 
field, (pix), for which 



Toi{x) = 7T{x)di(j){x) 



(46) 



with tt{x) = <p{x) the corresponding conjugate momentum. 

In the free scalar theory we can expand the energy- momentum tensor, and thus S{n), in 
terms of particle creation and annihilation operators. Employing the mode expansion of Toj 
in Eq. fl42p . we obtain for £^(n) in the free theory. 



lim R^ 

R—*oo 



dt 



d^p 



d^q 



(27r)32cup J (27r)32cu, 



-cUpU ■ q 



t p-«{'^p-'^q)*g«^"-(P-q) -|- h c 



(47) 
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We have chosen a normahzation for which [ap,ajj] = 2c<jp(27r)^(5^(p — q). We evaluate the 
integrals in spherical coordinates, performing the angular Op^q, (pp^g integrals in the stationary 
phase approximation. We find the stationary point 9p^q = On and 0p,g = that is, p, q are 
aligned with h. Then, 



1 /•°° q 

£{n) = lim — — — / dx~^ / dpp dqq — 



Sn«Le"'^"^""''^"^e-'^(P^-''^) (48) 



where we made use of the light-cone coordinates with respect to n = {l,n), = n-x = t—R, 
p"*" = n-p = ujp — p. The — > oo limit of the x~^ integral produces delta functions 
(27r)5(co'p — ujq) and {27T)6{ujp + ti;q). The latter delta function cannot be satisfied, so we are 
left with simply 

— — ^— u;p alapS^{n-p), (49) 

where we have normal-ordered the operators in Eq. ( I48l) and dropped the infinite constant 
term. The operator thus picks out the total energy of all scalar hadrons in the direction n. 
Inserting the result into Eq. fH3]) we reproduce the weighted cross section given in Eq. (1441) . 
completing the proof. 

There are corrections to this result from terms dropped in the stationary phase approxi- 
mation used in Eq. fHHl) of order 0{1/R), or, more precisely of order 0(l/(mo-R)), where mo 
is the smallest momentum scale that occurs naturally in the S-matrix elements. Although 
possibly small, it is still the inverse of some microscopic length scale. When multiplied by 
the macroscopic scale R, the result is a very large dimensionless number such that these 
corrections can safely be neglected. 

2. Dirac Fermions 

We can repeat the above derivation for Dirac fermions. The canonical energy-momentum 
tensor is 

T^,u = ii'lpidyip - g^^C. (50) 
Using the mode expansion of Toj in Eq. fj42|) gives for S{n) 

/"°° f d^p f d^q 
Sin) = lim R^ dt - — — / - — — n-a (51) 

^ ' «-°o Jo J (27r)32c^p J (27r)32a;q ^ ^ ^ 

X ^[6;,6^t?;"t(p)t;^(g)e-'("p-"'')*e'^" (P-'i) - a;,^a^u"^(p)M^(g)e'('"p-"'')*e-'^"-(P-'») 

r,s 

- 6;,a^t;^'^(p)u"(g)e~*('^p+'^'>)*e*^^-(P+'i) + a;,^6^^n't(p)t;^(g)e^('^p+^<')*e-^^"-(P+'i) 

Performing the angular integrals in the stationary phase approximation and taking the 
i? — i> oo limit of the t (or x^ = t — R) integral as in the previous section, we are left with 
the terms 

^{n) = ^^^^^3 t^PP— [aj,ftapn^i'""^(p^)^i'(p^) -&pft&J^'"^(p^)t^'(p^)]- (52) 
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The spinors satisfy the relations 

u''\pny{pn) = 2cJp(5"^ v^'^pnyipn) = 2up6''\ (53) 
and we obtain the final form of S{n), 

^^^^ = I (2^^'^^ Y.i^^'l + b5K)^\p - n), (54) 

after normal-ordering the 66^ term and dropping the infinite constant. In this form, S{n) 
picks out the energy of hadronic spin-1/2 Dirac fermions in direction n. 



B. Defining the thrust axis 

The factorization proof in Sec. [2] makes use of an operator t which gives the thrust axis 
t(X) when acting on a final state X. We begin by illustrating one explicit construction of 
t which we can use in full QCD. Then in SCET we argue that we can eliminate the use of 
the operator t entirely by identifying the thrust axis with the jet direction n in the two-jet 
current; an identification that is valid up to subleading power corrections in the event shapes 
we consider here. 



1. Thrust axis in QCD 

The thrust axis is the axis that is picked out by the maximum operation in the definition 
of thrust: ^ 

T = — max ^ |pi ■ t| . (55) 

i 

Each choice of axis t divides space into two hemispheres, A,B, so that the sum in Eq. 
splits into two corresponding pieces, 

T = ^ max {pA - Pb) ■ t , (56) 

where Pa,b are the total three-momentum in hemisphere A, B. By momentum conservation, 
Pa = -Ps, so 

T = max , (57) 

which is maximized by choosing 

t = r^, (58) 
IpaI 

the unit vector in the direction of the total three-momentum in hemisphere A. Thrust then 
simplifies to 

T = max ^^^^ . (59) 
A Q ^ ^ 
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Thus, to find the thrust, we need simply find the hemisphere with the largest total three- 
momentum, and the direction of this three-momentum is the thrust axis. To construct the 
operator which gives this axis, define the 4-vector of momentum fiow operators, 

POO 

V^{h) = (£{n),P{h)) = lim / dthiT^i{t,Rh) , (60) 

Jo 

extending the definition of the energy fiow operator S{h). Then the thrust axis operator is 

(61) 

where the quantity maximized by the max operator is the length of its three- vector argument, 
the integral is over the hemisphere A, and the operator J\f normalizes three- vectors, A/'(v) = 
v/ |v|. This construction manifests that the thrust axis of an event is just a function of 
particles' three-momenta and can be written in terms of the energy-momentum tensor in 
similar manner to the event shapes themselves. 



i = Af 



max / dQP(n) 



2. Thrust axis in SCET 

In SCET, we can replace t with the value of a "label" thrust axis whose determination 
is very simple. Namely, 11 is determined by the momentum labels in the two-jet operator 
Oniri2 ^-iid returns the thrust axis ti^{pi,p2) determined by these label momenta. The axis 
11 thus depends not on the final state X but only on the operator On^n2- For back-to-back 
jets, this "label" thrust axis is simply in the direction of the vector n in the operator Onn- 
In this section we argue that identifying this label thrust axis with the true thrust axis is 
valid up to power corrections that are subleading in the SCET expansion parameter A for 
angularities with a < 1. For a > 1 these corrections become leading order, and the thrust 
axis cannot be determined from label momenta alone. 

The true thrust axis t is that which maximizes the sum 

^ = ^El*-P^I' (62) 

where pj are the full three- momenta of all the particles in the event. The "label" thrust axis 
tL is defined to maximize 

^^ = ^El*L-Pi| , (63) 

where pj are the label momenta of all the particles in the event. The terms in Eq. (!63|) 
corresponding to each of the two jets may be grouped together, so that the individual label 
momenta in each group sum up to equal the total label momentum of each jet. This total 
label momentum is given by the label on the corresponding coUinear field in the operator 
Onin2) SO 11 depends not on the final state but only on this operator. In SCET, then, each 
of the operators e„, e^i, in Eq. (12^ depends on this axis 11, which in turn is completely 
determined by the choice of operator Cnin2- 

As we will show, the error induced in the event shape distribution by approximating the 
true thrust axis t by the label thrust axis 11 will be a subleading correction in the SCET 
expansion parameter A. The only condition is that the event shape in question sufficiently 
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suppresses the contribution of collinear particles close to the thrust axis, as is the case 
for angularity distributions with a < 1. To prove this, we will show that the two thrust 
axes are related by transformations of the light-cone vectors rii which leave the effective 
theory Lagrangian invariant, a property known as reparametrization invariance (RPI) j45| . 
They induce variations of the operators Onn and e appearing in da/de in Eq. (HM only at 
subleading order in A. 

The true and label thrust axes differ by a quantity 6t, where t = + (5t. The expression 
Eq. fl62|) that must be maximized to determine the true t is 

i^|t-(pi + k,)|, (64) 

where we expanded the momentum pj into its label and residual parts (for soft particles 
pi = 0). This expression is identical, through terms of order A, to the expression that must 
be maximized to determine the label thrust axis 11, Eq. fl63l) . Thus t = 11 through order A, 
and we conclude that St is of order A^. Furthermore, since t and are both unit vectors, 
we have 

1 = \tf = |tL + 6tf = \ti^f + 2tL ■ 6t + \6tf = 1 + 2tL ■ 5t + O(A^) . (65) 

Thus 11 ■ 5t = + (9(A^); that is, 6t is orthogonal to 11 up to terms of order A^. For a 
two-jet operator Onn, the label thrust axis 11 is just n, and 6t is purely transverse to n: 

Due to RPI, each of the collinear sectors in SCET is invariant under changes in the light- 
cone vectors rzj and nj that leave the conditions = nf = and rii-ni = 2 unchanged. Here, 
we are considering an arbitrary number of collinear sectors each labeled by i, and associated 
with a corresponding light-cone direction rii, where = (1, nj),nf = (1, — rij) with = 1. 
Note that for only two label directions, momentum conservation fixes = ni^ = and 

= = ^'^^ however, for the analysis below we will continue to distinguish i = 1,2 until 
the very end. 



Under RPI type-I and type-II transformations [45| in the ith. collinear sector the four- 



vectors Hi, rii are shifted by a transverse pieces A^,ej 

(I) ni^ni + (II) ni^ni + , (66) 

where is allowed to be 0{X) or smaller, and ej- is allowed to be 0{1) or smaller. Below 
we will choose both A^,e^ to be only (9(A^). Under these transformations label momenta 
= Ui-pn^ /2 +p^^ in the zth sector transform as 

(I)p/^^p'/^=pM_lA^^.p^nf (67a) 

{ll)P''^P"' = P'-\ei-p^nt. (67b) 

Since the label thrust axis tL depends on label momenta in each collinear sector, tL 
also transforms under reparameterizations of nj,nj in any sector: 11 — > tL' = tL + ^tj^. 
Thus we can bring 11 to coincide with the full thrust axis t by performing a suitable set 
of RPI transformations on nj,nj. For the case of two jets, we will use an RPI type-I and 
a type-II transformation of each of ni and n2 to bring 11 to equal t. To find the correct 
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transformations, we require that the full thrust axis t maximize the sum in Eq. (163|) after 
the label momenta are transformed as in Eq. (!67|) . After these transformations, the sum 
takes the form 



E 



X ■ p' 



+E 

ie2 



X 



Pi 



1 

2' 



2 
1 



(68) 



where the sum is divided into the rii and n2 sectors, and x is the variable vector that 
must be chosen to maximize the sum. Let us use that rii = —112 = n, and consider RPI 



transformations such that A]^ 
Eq. (EH]) as 



E 



Aj_ and ef 



(x-n)(Ax-£±) 



Pi 



Then we can write 
(69) 



But this is the same form as the sum giving the label thrust axis in Eq. (!63l) . so we know 
that Eq. (l69l) is maximized by choosing x such that the vector in brackets is just 11 itself. 



-(x-n)(A^ -£j 



(70) 



Using our choice of power counting, A±,e± ~ 0{X^), we can solve for x order-by-order in 
A. The solution starts as x = 11 + C(A^), and we can use • n = 1, to obtain 



X = tL + -(Ax-£x). (71) 

This new solution for the label thrust axis coincides with the true thrust axis if we choose 

A± = -e± = 6t . (72) 
In terms of four-vectors, Eq. ( l72l) corresponds to the set of RPI transformations, 



ni 
n2 



;i,ni + <5t) 
;i,n2-5t) 



ni (1, —111 — (^t) 
n2 (1, -n2 + 5t) . 



(73) 



The label thrust axis corresponding to the operator On^n2 is thereby brought to coincide 
with the true thrust axis. These leave the SCET Lagrangian invariant. The operator Onn 
is invariant up to corrections subleading in A. 

Finally, we need to estimate the size of the variation in the operator e given in Eq. f|T3l) . 
but with the full thrust axis replaced by the label thrust axis, which is just the direction 
n appearing in Onn- Consider again the example of angularities e = r^, for which fa{fl) = 
g-|';l(i-a)^ following As in Eq. ( l20l) . the contribution of an individual particle to Ta{X) is 



iPrI 

Q 



-H(l-a) 



IPtI 

Q 



mm 



Et_ 
E± 



l-g 
2 



(74) 



We recall from Sec. 12 DI that for usoft particles pt ~ QA^ and the ratio in parentheses i? ~ 1, 
while for coUinear particles px ~ QA and ~ A^. The shift in the thrust axis identified 
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above induces for usoft particles 5pT ~ QA^ and 5R ~ A^, while for collinear particles we 
find 5pT ~ QA^ and 5R ~ A^. This changes the total contribution of a usoft particle to e{X) 
by 5es ~ A^, and of a collinear particle by 5ec ~ A^~". 

As long as a < 1, Se^ is smaller than A^. If, as in this paper, we are interested in 
calculating the event shape accurately only to order A^, then we may neglect the shifts 
and Scc- Beginning with a = 1, for which ti = B (broadening), the power corrections 
induced by the shift in the thrust axis cannot be neglected, since Scc ~ A^~" becomes as 
large as the terms we considered in this paper. The necessity of accounting for this "recoil" 
of the thrust axis against jisoft-scale momenta was demonstrated in full QCD for B in 37 
and for r„ with a > 1 in 



4. CONCLUSIONS 

We have proved the factorization of two-jet event shape distributions in soft collinear 
effective theory without assuming the factorization of hadronic final states into separate 
colored collinear and soft sectors as in previous discussions. To do this we expressed the 
weight of each final state 6{e — e{X)) in the differential cross section da/de as an operator 
built out of energy flow and thrust axis operators acting on the final state, allowing us to sum 
over a complete set of hadronic states before factorizing soft and collinear matrix elements. 
These results are valid up to 0{X) corrections to the decoupling of usoft and collinear degrees 
of freedom in SCET, 0{1/R) corrections to the relation between the energy flow operator 
and energy-momentum tensor, and (9(A^) corrections due to the difference between the 
thrust axis and the collinear jet direction n, where the power b depends on the event shape 
in question. Similar methods should also be useful in studying the factorization of other jet 
observables in both leptonic and hadronic collisions. 
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